Critical Currents in Quasiperiodic Pinning Arrays: 
Chains and Penrose Lattices 
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We study the critical depinning current J c versus the applied magnetic flux "!>, for quasiperiodic 
(QP) chains and 2D arrays of pinning centers placed on the nodes of a five-fold Penrose lattice. In 
QP chains, the peaks in J c (<&) are determined by a sequence of harmonics of the long and short 
segments of the chain. The critical current J c (&) has a remarkable self-similarity. In 2D QP pinning 
arrays, we predict analytically and numerically the main features of J c ($), and demonstrate that 
the Penrose lattice of pinning sites provides an enormous enhancement of J c ($), even compared 
to triangular and random pinning site arrays. This huge increase in Jc($) could be useful for 
applications. 

PACS numbers: 74.25.Qt 



Recent progress in the fabrication of nanostructures 
has provided a wide variety of well-controlled vortex- 
confinement topologies, including different regular pin- 
ning arrays. A main fundamental question in this field 
was how to drastically increase vortex pinning, and thus 
the critical current J c , using artificially-produced peri- 
odic Arrays of Pinning Sites (APS) UH, HSU Tne 
increase and, more generally, control of the critical cur- 
rent J c in superconductors by its patterning (perforation) 
can be of practical importance for applications in micro- 
electronic devices. 

A peak in the critical current J c ($), for a given value 
of the magnetifc flux per unit cell, say $i, can be engi- 
neered using a superconducting sample with a periodic 
APS with a matching field Hi = &i/A (where A is the 
area of the pinning cell), corresponding to one trapped 
vortex per pinning site. However, this peak in J C ( ( E > ), 
while useful to obtain, decreases very quickly for fluxes 
away from <I>i. Thus, the desired peak in J C ( < I ) ) is too nar- 
row and not very robust against changes in $. It would 
be greatly desirable to have samples with APS with many 
periods (ideally infinite). This multiple-period APS sam- 
ple would provide either very many peaks or an extremely 
broad peak in J C (<E>), as opposed to just one (narrow) 
main peak (and its harmonics). We achieve this goal [a 
very broad J c (<&)] here by studying samples with many 
built-in periods. 

Here, we study vortex pinning by quasiperiodic (QP) 
chains and by 2D APS located on the nodes of QP lat- 
tices (e.g., a five- fold Penrose lattice) @. We show that 
the use of the 2D QP (Penrose) lattice of pinning sites re- 
sults in a remarkable enhancement of J c (<!>), as compared 
to other APS, including triangular and random APS. In 
contrast to superconducting networks, for which only the 
areas of the network plaquettes play a role , for vortex 
pinning by QP pinning arrays, the specific geometry of 
the elements which form the QP lattice and their arrange- 
ment (and not just the areas) are important, making the 
problem far more complicated. 

Simulation. — We model a three-dimensional (3D) 



slab, infinitely long in the z-direction, by a 2D (in the 
xy-plane) simulation cell with periodic boundary con- 
ditions. We perform simulated annealing simulations 
by numerically integrating the overdamped equations 
of motion (see, e.g., Ref. @, 0): rrvj = f, = 

fvv + fvp + f T + f d Here f . ig the tQtal force per umt 

length acting on vortex i, f™ and f" p are the forces 
due to vortex-vortex and vortex-pin interactions, re- 
spectively, fj T is the thermal stochastic force, and ff 
is the driving force; rj is the viscosity, which is set to 
unity. The force due to the vortex-vortex interaction is 
fT = Ef" fo Ti - rj | /A)fy, where N v is the 
number of vortices, K\ is a modified Bessel function, A is 
the penetration depth, f,j = (r— rj)/ | r,*— rj |, and fo = 
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where N p is the number of pinning sites, f p (expressed 
in fo) is the maximum pinning force of each short-range 

parabolic potential well located at r). , r p is the range of 
the pinning potential, is the Heaviside step function, 
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All the lengths (fields) 



are expressed in units of A («f>o/A 2 ). The ground state of 
a system of moving vortices is obtained by simulating the 
field-cooled experiments (e.g., [13 )• For deep short-range 
(5-hke) potential wells, the energy required to depin vor- 
tices is proportional to the number of pinned vortices, 
Ny^ . Therefore, in this approximation, we can define 
the critical current as follows: j c ( ( f ) ) = joNi P \&)/N v (&), 
where jo is a constant, and study the dimensionless value 
Jc = jc/ja- We use narrow potential wells as pinning 



sites, with 



0.04A to 0.1A. We have also performed 



dynamical simulations of J c using a threshold criterion, 
i.e., J c is obtained as the minimum current J oc ff 
which depins the vortices. The results obtained using 
these two criteria are essentially equal [ll|, and here we 
use the "static" criterion defined above. 

ID quasicrystal. — A QP chain 6] can be con- 
structed by iteratively applying the Fibonacci rule (L — > 
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FIG. 1: (Color online) (a) The critical depinning current J c , 
versus the number of vortices, N v ~ $, for ID QP chains, 
N p = 21 (red bottom line), N p = 34 (second from the bottom 
blue line), N p = 55 (green line), and N p — 89 (dark blue 
top line), for 7 = as/a.L = 1/t. Here we use: f p /fo = 1.0 
and r p = 0.1 A. Independently of the length of the chain, the 
peaks include the sequence of successive Fibonacci numbers 
and their sub-harmonics, (b) The function J c ($/$i) for the 
same ID chains (using same colors). The curves for different 
chains display the same set of peaks, namely, at = 1 

(first matching field) and <J>/<I>i = 0.5, as well as at the golden- 
mean- related values: $/$i = r, r/2, (r + l)/2 = r 2 /2, (r 2 + 
t)/2 = t 3 /2, t 2 = T + 1, r 2 + 1. This behavior demonstrates 
the self-similarity of J c ($). 

LS, S — > L), which generates an infinite sequence of 
two line segments, long L and short S. For an infinite 
QP sequence |(| , the ratio of the numbers of long to short 
segments is the golden mean r = (1 + VE)/2. Let, e.g., 
as = 1 and = r, where as and at are the lengths 
of the short and long segments, respectively. Then the 
position of the nth point where a new segment, either L 
or S, begins is determined by: x n = n + [n/r] /r, where 
[x] denotes the integer part of x. To study the critical 
depinning current J c in QP pinning chains, we place pin- 
ning sites to the points where the L or S elements of the 
QP sequence link to each other. 

The results of calculating J C (N V ) for chains of different 



length and the same 7 = 1/r are shown in Fig. la. The 
plot clearly shows that, for sufficiently long chains, the 
positions of the main peaks in J c , to a significant extent, 
do not depend on the length of the chain. The peaks 
form a Fibonacci sequence: N v = 13, 21, 34, 55, 89, 144, 
and other "harmonics": N v = 17, 27.5 (= 55/2), 44.5 
(= 89/2), etc. Of course, longer chains allow to better 
reveal peaks for larger Fibonacci numbers. In Fig. lb, the 
same curves are rescaled, normalized by the numbers of 
pins in each chain. The rescaled J c curves reproduce each 
other and have many pronounced peaks for golden-mean- 
related values of <£>/$i ($1 is the flux corresponding to 
the first matching field, Hi, when N v = N p ), as shown in 
Fig. lb. Therefore, the same peaks of J c (&), for different 
chains, are revealed before and after rescaling because of 
the self-similarity of J c (®)- The self-similarity of J c (&) 
has also been studied in reciprocal A:-space and will be 
presented elsewhere [Tl| . 

Penrose lattice. — Consider now a 2D QP APS, 
namely, an APS located at the nodes of a five-fold Pen- 
rose lattice. This lattice is a 2D QP structure, or qua- 
sicrystal, also referred to as Penrose tiling 0. These 
structures possess a perfect local rotational (five- or ten- 
fold) symmetry, but do not have translational long-range 
order. The unusual self-similar difraction pattern of a 
Penrose lattice exhibits a dense set of "Bragg" peaks be- 
cause the lattice contains an infinite number of periods 
in it 0. It is precisely this unusual property that is re- 
sponsible for the striking J c (<E>)'s obtained here. 

The structure of a five-fold Penrose lattice is presented 
in Fig. 2a. As an illustration only, a small five-fold sym- 
metric fragment with 46 points is shown. According to 
specific rules, the points are connected by lines in order 
to display the structure of the Penrose lattice. The el- 
emental building blocks are rhombuses with equal sides 
a and angles which are multiples of 8 — 36°. There are 
two kinds of rhombuses: (i) those having angles 29 and 
36* (so-called "thick"; shown empty in Fig. 2a), and (ii) 
rhombuses with angles 9 and 46* (so-called "thin" ; filled 
in orange in Fig. 2a). 

Let us analyze whether any specific matching effects 
can exist between the Penrose pinning lattice and the 
interacting vortices, which affect the magnetic-field de- 
pendence of the critical depinning current J c (&) (Fig. 2). 
Quasicrystalline patterns are intrinsically incommensu- 
rate with the flux lattice for any value of the magnetic 
field |7J, therefore, in contrast to periodic (e.g., triangu- 
lar or square) pinning arrays, one might a priori assume 
a lack of sharp peaks in J c ($) for QP APS. However, the 
existence of many periods in the Penrose lattice can lead 
to a hierarchy of matching effects for certain values of 
the applied magnetic field, resulting in strikingly-broad 
shapes for J c {^)- These could be valuable for applica- 
tions demanding unusually broad J c ($)'s. 

To match the vortex lattice on an entire QP APS, the 
specific geometry of the elements which form the QP lat- 
tice is important, as well as their arrangement. While 
the sides, a, of the rhombuses are equal, the distances 
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FIG. 2: (Color online) (a) An example of a five-fold Pen- 
rose lattice, consisting of "thick" and "thin" rhombuses, (b 
to d) The location of vortices (green dots) and the Penrose- 
lattice APS (red open circles connected by black lines) for: 
(b) $ = $vacancy/thin = "^v/t = 0.757 $1, vortices occupy 
all the pinning sites except those in one of the two vertices 
of each "thin" rhombus; (c) $ = $1, N v coincides with N p ; 
(d) $ = $ ln terstitial/thick = *i/T = 1.482 $1, vortices occupy 
both the pinning sites and the interstitial positions inside each 
"thick" rhombus. Here f p /f = 2.0, r p = 0.1A. (e) The 
J c ($ ~ N v ) for Penrose-lattice arrays with N p = 46, 141, 301. 
The data points [b], [c], [d] refer to snapshots in 2b, 2c, 2d, re- 
spectively. The peak [c] at $i is suppressed for weaker pinning 
(magenta dotted J c curve) . Eventually, all the main peaks dis- 
appear for sufficiently weak pins (black dashed curve). The 
inset shows the dimensionless difference, /diff, of the pinning 
and the elastic energies versus the pinning-to-interaction en- 
ergy ratio, for the broad J c peak at $ v / t (red dashed line) 
and for 3>i (red solid line). Only / d iff > gives stable peaks 
in J c . (f) The J c (<&) for a 301-sites Penrose-lattice (red solid 
line), (recalculated for flux only on the Penrose area, Ap), 
triangular (black solid line) and random (green open circles 
and solid line) pinning arrays. 



between the nodes are not equal. The lengths of the diag- 
onals of the rhombuses (Fig. 2a) are: 1.176a, ra w 1.618a 
(for the thick rhombus), (r — l)a = a/r w 0.618a, and 
1.902a (for the thin rhombus). Based on these distances, 
we can predict matching effects (and corresponding fea- 
tures of the function J c ($)) for the Penrose-lattice APS. 

First, there is a "first matching field" (we denote the 
corresponding flux as $i) when each pinning site is oc- 
cupied by a vortex (Fig. 2c). Although the sides of all 
the rhombuses are equal to each other, nevertheless this 
matching effect is not expected to be accompanied by a 
sharp peak. Instead, it is a broad maximum (peak [c] in 
Fig. 2e) involving three kinds of local "commensurabil- 
ity" effects of the flux lattice: with the rhombus side a; 
with the short diagonal of a thick rhombus, 1.176a, which 
is close to a; and with the short diagonal of a thin rhom- 
bus, which is a/r ss 0.618a. For the overall square cell 
used, some of the vortices are outside the Penrose sample; 
these mimic the applied magnetic field and determine the 
average vortex density in the entire cell. Because of the 
additional vortices outside the Penrose APS, our com- 
puted J c (&) is reduced by a factor r\ = Ap/A w 0.575 
(see Fig. 2e), where Ap and A are the areas of the Pen- 
rose lattice and of the cell. 

Another matching (Fig. 2b) is related with the filling 
of all the pinning sites on the vertices of the thick rhom- 
buses and only three out of four of the pinning sites on 
the vertices of thin rhombuses, i.e., one of the pinning 
sites on the vertices of the thin rhombuses is empty. For 
this value of the flux, matching conditions are fulfilled for 
two close distances, a (the side of a rhombus) and 1.176a 
(the short diagonal of a thick rhombus), but are not ful- 
filled for the short diagonal, a/r, of the thin rhombus. 
Therefore, this 2D QP feature is related to t, although 
not in such a direct way as in the case of a ID QP pin- 
ning array. This 2D QP matching results in a very wide 
maximum (arrow [b] in Fig. 2e) of the function J c (&). 
The position of this broad maximum (denoted here by 
$vacancy/thin = *v/t = 0.757 $i) could be found as fol- 
lows. The ratio of the numbers of thick and thin rhom- 
buses is determined by the Fibonacci numbers and in the 
limit of large pinning arrays, N p — ► oo, this ratio tends to 
t. The number of unoccupied pinning sites is governed 
by the number of thin rhombuses. However, some of the 
thin rhombuses are separated from other thin rhombuses 
by a single thick one (single thin rhombuses; see Fig. 2a), 
while some of thin rhombuses have common sides with 
each other (orange arrow-shaped double thin rhombuses 
in Fig. 2a). Therefore, the number of vacancies (i.e., 
unoccupied pins) is then the number of single thin rhom- 
buses N* h plus one half of the number N^ h of "double" 
thin rhombuses, A^ n ($ v/t ) = A^ h + A^/2, where A^ n is 
the number of unoccupied pinning sites at $ = $ v / t - 

For higher vortex densities ($ = ^interstitial/thick = 
$i/T = 1-482 <I>i) a single interstitial vortex is inside each 
thick rhombus (see Fig. 2d). These interstitial vortices 
can easily move; thus J c has no peak at $i/T- The posi- 
tion of this feature is determined by the number of vor- 
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tices at $ = $1, which is N v ($) = N p , plus the number tic energy with respect to the pinning energy; the latter 



of thick rhombuses, N^ ick = N th /r. 

In order to better understand the structure of J c (&) 
for the Penrose pinning lattice, we compare the elastic 
E c \ and pinning E p i n energies of the vortex lattice at 
Hi and at (the lower field) H v u, corresponding to the 
two maxima of J c (Figs. 2e,f). Vortices can be pinned 
if the gain E pm = Up- m n pm of the pinnin g e nergy is 

related 



larger than the increase of the elastic energy 



to local compressions: E c \ = Cn[(c 



b)/a c 



The 



shear elastic energy (ex C 66 ) provides the same qualita- 



tive result Here, £/, 



pin 



fp r pj 



kpin 



is the density 



of pinning centers, 0(H < Hi) = H/Hi = B/($ n p i n ), 
and 0(H > Hi) = 1 is the fraction of occupied pin- 
ning sites (0 = 1 for H = Hi, and — 0.757 for 

1/2 

H = H v / t ), a oq = (2/v / 3/3n p in) is the equilibrium 
distance between vortices in the triangular lattice, b 
is the minimum distance between vortices in the dis- 
torted pinned vortex lattice (b = a/r for H = Hi and 
b = a for H = H v/t ), and C u = S 2 /[4tt(1 + X 2 k 2 )} 
is the compressibility modulus for short-range deforma- 
tions ^2] with characteristic spatial scale k w (npin) . 
The dimensionless difference of the pinning and elas- 
tic energies is E pin - E e \ = 0f di s 7i pin $g/(47rA 2 ), where 

/diff = 47rA 2 C/ pin /$ 2 1 - [l - b (/?V3n pin /2) 1/2 j ' . Near 

matching fields, J c has a peak when /diff > (and no 
peak when /diff < 0). Since only two matching fields 
provide /diff > 0, then our analysis explains the two- 
peak structure observed in J c shown in Figs. 2e,f. For 
instance, for the main matching fields: /diff(3?v/t) ~ 
0.0056, /diff(*i) » 0.0058, and /diff($i /T ) « -0.09. 
Note that for weaker pinning, the two-peak structure 
gradually turns into one very broad peak, and eventu- 
ally zero peaks for weak enough pinning (see Fig. 2e). 
The J c peaks corresponding to higher matching fields are 
strongly suppressed because of the fast increase (oc B 2 ) 
of the compressibility modulus Cu and, thus, the elas- 



cannot exceed the maximum value U p \ n n pin . The sub- 
harmonic peaks of J c , which could occur for lower fields 
H < .ffv/t! are a l so suppressed due to the increase of 
Cu associated with the growing spatial scales 1/k of the 
deformations. 

For comparison, we show the J c (&) for the Penrose- 
lattice (itself, i.e. calculated only for Ap), triangular 
and random pinning arrays (Fig. 2f). The latter is an 
average over five realizations of disorder. Notice that the 
QP lattice leads to a very broad and potentially useful 
enhancement of the critical current J c ($), even compared 
to the triangular or random APS. The remarkably broad 
maximum in J c ($) is due to the fact that the Penrose 
lattice has many (infinite, in the thermodynamic limit) 
periodicities built in it [6j. In principle, each one of these 
periods provides a peak in J C (<1>). In practice, like in 
quasicrystalline difraction patterns, only few peaks are 
strong. This is also consistent with our study. Further- 
more, the pinning parameters can be adjusted by using 
as pinning centers either antidots "drilled" in the film 
0,0, or blind antidots of different depths and radii. 
Thus, our results could be observed experimentally. 

Conclusions. — The critical depinning current J c {^) 
was studied in ID QP chains and in 2D QP arrays (the 
five-fold Penrose lattice) of pinning sites. A hierarchical 
and self-similar J c (&) was obtained. We physically ana- 
lyzed all the main features of J C ( < I > ). Our analysis shows 
that the QP lattice provides an unusually broad critical 
current J C ( < I J ), that could be useful for practical appli- 
cations demanding high J c 's over a wide range of fields. 
Our proposal can be easily extended, mutatis mutandi, 
to other related systems, including colloidal suspensions 
interacting with pinning traps provided by arrays of op- 
tical tweezers [l3j . 
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